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SUMMARY 

The pump-source intensity distribution inside a cylindrical dye cell has appli- 
cation to several aspects of dye-laser operation: (1 ) internal power and energy 

distribution, (2) lasing and coupling efficiency, (3) thermal and optical perfor- 
mance, and (4) pump-source selection. The main impetus for the work presented here 
is pump-source selection. Specifically, the internal intensity and absorbed-powe r 
distribution of a simplified hypothetical dye laser of cylindrical geometry is 
calculated. Total absorbed power is also calculated and compared with laboratory 
measurements of lasing-threshold energy deposition in a dye cell to determine the 
suitability of solar radiation as a pump source or, alternatively, what modifica- 
tions, if any, are necessary to the hypothetical system for solar pumping. 


INTRODUCTION 

NASA research on laser power transmission in space began in 1972 and involved 
several types of lasers (ref. 1). In the late 1970's the program focused on the use 
of solar-pumped lasers. Several prospective active media for solar pumping have been 
identified and investigated in the program. One of these is dye. Dyes offer com- 
paratively good solar spectral absorption, lasing efficiency, excellent tunability, 
and low cost. However, very large light intensities are required for their opera- 
tion. A question of prime importance then is whether focused sunlight can pump dye 
lasers. The theory reported here, although available for broader application, is 
used with experimental data to provide information relevant to this question. 
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SYMBOLS 

area element 

area element of dye-tube surface 
solar-concentration factor 
solar spectral power distribution 
intensity from a particular direction 

2 

solar constant, W/mm 

absorbed intensity from particular direction 
defined in equation (6) 

fluorescent frequency conversion factor 
geometric absorption factor 
length factor 
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dye-molecule factor 


' ^2 


2 

intensity inside cylinder, W/mm 

2 

intensity on area element, W/mm -sr 
concentration 

power emitted into interior, W 

power absorbed from all directions 

power absorbed from particular direction 

radius from center of cylinder to dA element 

radius of cylinder, mm 

distance between dA g and dA, mm 

Fresnel transmission factor 

volume 

dye -absorption coefficient, mm -1 

wave length -dependent absorption coefficient 

angle of solar subtense, radians 

angle between Rq and R, radians 

wavelength 

lower and upper bounds of wavelength, respectively 
spectral molecular-absorption cross section 
angle between Rq and , radians 


MODEL THEORY 

A long cylindrical geometry for the dye cell is postulated for the model. The 
cell is surrounded by an optical concentrator which collects solar radiation in space 
through a large aperture and directs central rays from the Sun at the cell axis per- 
pendicular to it. The concentrator surface is distant enough from the cell that from 
any point on the surface of the concentrator, the cell subtends a smaller angle (nor- 
mal to the cell axis) than that of the solar disc. Such concentrators are possible. 
Further consideration of the concentrator is not germane to the model presented 
here. Thus, the dye tube is bathed with uniform concentrated sunlight over its 
length. At any point on the surface of the dye tube, light is incident through 180° 
in any plane normal to the tube axis. In all planes containing the tube axis, light 
is incident at any point on the tube surface at angles of 0° through 0.26° (half the 
angle of solar subtense y). With this idealized but feasible geometry, intensity 
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within the dye cell will have only radial dependence. With term definitions given in 
figure 1, the equation for intensity within the dye cell can be developed as follows: 

With Iq incident on the area element of the dye-tube surface dA g , the 
differential power emitted into the interior is 


2 dA , „ . 

d P = L dA — - cos 4> O ) 

0 s r 2 

The differential intensity at the area element dA in the dye is 

-ctR dA 

dl = I Q Te — 2^ COS ^ (2) 
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T Fresnel transmission factor 

a dye -absorption coefficient , mm' 
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Expressing I Q in terms of the solar constant E g and solar-concentration factor C 
(ref. 2) gives 
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( 6 ) 


The following approximations have been included in the derivation: 

(1) Scattered light is negligible. 

(2) Dimensions of the concentrating optics are much greater than the 

tube radius. 

Now consider a particular radial position in the dye. Rays of light converge 
there from 360° in a plane normal to the axis. The intensity of any particular ray 
(denoted by a prime) undergoes a differential absorption in a distance dR^ given by 

dE^ = a E* (0) dR 1 (7) 

The power absorbed from a particular direction is 

A dE' = a E* (0) A dR (8) 

a 1 

a E’ (0) (9) 


dP* = 
a 


dP' 
a 

dV 


Therefore, the differential power absorbed per unit volume from all directions is 
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dV 


a / E' (0) d0 = aE C F(a,R) 
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Since attenuation is wavelength dependent, 


dP 
a 

dV 



a(\) F(a,R) 


dE 




a(\) E(X) F(a,R) dX 
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where E(\) is the solar spectral power distribution. 


4 



CALCULATIONS 


In equation (6), F(a,R) was evaluated at 21 radial positions (from centerline 
to periphery, 2 mm) and with 11 dye concentrations (absorption coefficient, 

0 to 2 mm” 1 ). (It should be noted that F(a,R) may be normalized to F(ccRq,R/Rq) 
for application to other laser geometries, ) The calculations include assumptions 
that: (1) the index of refraction of the dye solution is 1.32, (2) the index of 

refraction of the glass tube enclosing the dye is 1.5, and (3) reflection at the 
interface of the glass and solution is negligible. Each data point represents the 
sum of incremental intensities from every direction in a plane normal to the dye-tube 
axis. The increment used for the summation was 0.5°, which is approximately the 
angle subtended by the Sun. Actual data values are presented in table I in groups 
of rows and columns. Each group corresponds to a particular dye -absorption coef- 
ficient a; from top to bottom the coefficients are a = 0, 0.2, 0.4, 0.6, ..., 

2.0 mm"" 1 , respectively. Within each group the value in the upper left corner is the 

normalized intensity at the dye-tube centerline (R/Rq = 0). Successively larger 
values of R/Rq (0*05, 0.10, 0.15, ...) are read from left to right, top to bottom, 

in each group. 

The data of table I are plotted in figure 2. The top curve corresponds to 
a = 0 (no dye present). Successively lower curves correspond to a = 0. 2 through 

2.0 mm" 1 . 

Spectral values of a(\) and E(\) were taken from references 3 and 4, respec- 
tively, and used with values of F ( a, R) interpolated from table I to calculate 
absorbed-power density dP a /dV at 21 radial positions. Figure 3 shews dP a /dV 
plotted against radial position R/Rq. From top to bottom, the curves correspond 
to 4n, 2n, n, and n/2, respectively, where the concentration n is 

1.51 x 10 13 molecules per cubic millimeter, the dye concentration of the experiment 
discussed in the appendix. (Here, a(\) varies with n according to a(\) = n a(\), 
where a(\) is the spectral molecular-absorption cross section.) 

The absorbed-power density dP a /dV integrated over the volume of the tube gives 
the total power absorbed by the dye. The integration was performed numerically for 
all curves of figure 3 and was evaluated at various dye concentrations as follows: 


Concentration 

Power, W 

n/2 

0.02623C 

n 

. 038C 

2n 

• 052C 

4n 

. 066C 


DISCUSSION 

The curves of F { a, R ) in figure 2 show expected trends at the centerline and 
boundary of the dye. At the centerline (R/Rq = 0), intensity decreases exponentially 
according to 


F( a, 0) 



( 12 ) 
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as dye concentration n increases. This is caused, of course, by higher concen- 
trations of dye absorbing more light in the outer portions of the dye before the 
light can reach the center of the tube. The integral form of the normalized 
intensity (see eq. (6)) reduces exactly to the exponential form in equation (12) 
at the centerline where = Rq, T = 0.96, and <J) = 0. At the outer boundary 

(R/R Q = 1) the integral equation does not reduce to a simple form, but there is a 
more rapid attenuation toward an intensity limit which is determined by the amount of 
light that can be gathered through the tube wall. (In the limit as a becomes very 
large, negligible light would come from the dye side of the outer boundary.) Calcu- 
lations of the intensities at the outer boundary were not done because of time limi- 
tations associated with the coirputer program. Rather, those intensities were 
extrapolated from the nearest three data points of each curve. 

At radial positions between the two extremes, the most notable feature is the 
abrupt intensity decrease near R/R Q = 0*7. That feature is caused by Fresnel 
reflection at the surface of the glass tube. At these radial distances a range of 
0 values exists for which the corresponding refraction angles are impossible, 
meaning that no light can reach some of the larger radial positions from some points 
on the surface of the dye tube. The greater the radial position, the greater these 
"blind" and "nearly blind" areas are on the tube surface. 

The data of table I are sufficient to interpolate F(a,R) for a values only 
as large as 2. This limits the number of curves in figure 3. If more curves had 
been calculated, their character would vary approximately as indicated in figure 4; 
that is, values near the centerline would begin to decrease, and values at the 
boundary would continue to increase. 

None of the curves of figure 3 have zero slope over extended lengths. This 
means that there is no absolutely uniform energy deposition over large areas. 

(Absolute uniform deposition is achieved only for the trivial case where there is no 
energy deposition at all.) This indicates that a radial thermal gradient will always 
exist in such a laser volume with such irradiation. Gradient magnitudes will vary 
during illumination until an equilibrium is established and will cause variations in 
the refractive index of the fluid which, in turn, cause lens like effects in the 
optical cavity that can destroy optical feedback. Longer cavities are more sensitive 
to the gradients if other laser parameters remain unchanged. The growth rate of the 
gradients can be lessened by smaller dye concentrations, in which case a longer 
length is needed to compensate for loss of optical gain. Steady -state gradients 
could be minimized by optically designing solar radiation to be focused more toward 
the center of the tube. If, simultaneously, dye concentration were increased, total 
power deposition would be increased. 

The absorbed-power density dP a /dV integrated over the cylindrical volume of 
the dye was evaluated previously. This information may be combined with data from 
our laboratory dye -laser system (see the appendix) to provide an estimate of the 
solar-concentration factor required to operate a solar-pumped dye laser (like our 
laboratory system) in space. 

If this same laboratory laser were illuminated in space as our model theory 
describes, the total power absorption as calculated from theory approximates the 
absorbed power that was measured in the laboratory (26 700 W) in order to reach 
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lasing threshold* From this relationship an estimate of the solar concentration 
required to reach lasing threshold is 


C 


26 700 
0.038 


7x10 
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This is about 35 times the maximum geometrical solar concentration attainable in 
practice. (The Sun does not approximate a point source well enough to achieve 
practical concentrations greater than 20 000. ) Therefore, a feasible model must 
incorporate irnprovements that compensate for the factor 35 while maintaining a solar 
concentration of 20 000. Improvements can take several forms, some of which are 
stated as follows: (1) absorption by the dye mixture can be made more efficient 

both geometrically and spectrally, (2) fluorescent frequency conversion may be used, 
(3) the rhodamine-dye molecule may be improved, and (4) laser-gain length can be 
increased. Each of these possible improvements will be addressed briefly. 

Recall the absorbed powers calculated previously. These data may be extrapo- 
lated to very high dye densities to show that absorbed power may be increased by a 
factor of 1.7 by increasing dye concentration. Conversely, this same factor could 
be gained geometrically by multipass absorption in the same dye concentration rather 
than by the single-pass absorption of the model. This factor is labeled the 
"geometric absorption factor" F • 

Outside the absorption band of all dyes, fluorescent frequency conversion could 
also be used to translate incoming radiation into absorption bands (refs. 5 and 6). 
(Fluorescent frequency conversion may be viewed as nongeome trie solar concentration. ) 

From a theoretical treatment of dye-laser threshold, reference 7 shows that the 
lasing threshold of the rhodamine 6G molecule cannot be reduced by more than a factor 
of 4. This reference also shows that self -tuned critical inversion is inversely 
proportional to the three-fourths power of laser length (if triplet population is 
kept small and/or cavity losses are sufficiently large). 

The product of the aforementioned factors must be greater than 35. That is. 


F F F F > 35 
g f m x 

where 

F^ fluorescent frequency conversion factor 

F dye-molecule factor 

m 

F^ length factor 

If we estimate the factors as follows: 

F = 1 .5 F = 3 F = 1 

g f m 
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and 


F 




> 


35 

4.5 


7.77 


then would correspond to a length change by a factor of 15.41. That is, the 

active length of the model laser system must be 4.72 m to achieve a lasing 
threshold. It is practically impossible to build a solar concentrator that will 
provide the radiation geometry of this lengthened model at a solar concentration of 
20 000. Therefore, it is not possible to build a solar-puirped dye laser of this 
particular lengthened geometry in practice. However, another geometry that provides 
the increased length (with essentially the same laser cavity and absorbed -power 
density) might be possible, assuming, of course, that threshold-length scaling 
actually follows theoretical predictions at the required length. 


CONCLUDING REMARKS 

The intensity and power distribution inside a cylindrical dye cell have been 
calculated and presented in normalized form for application to a model system and to 
other similar systems. Also, the relevant equations were derived and presented for 
use on systems with major geometrical differences (e.g., nonuniform illumination). 
These equations and calculations serve as a basis for investigations into other areas 
of dye-laser research, such as thermal, thermo-optical, modal, and optical analyses. 

For the particular space-based model examined here, analysis shows that radial 
gradients in absorbed-power deposition exist and increase with dye concentration. 
(Radial gradients cause lenslike effects deleterious to lasing.) Total-power deposi- 
tion also increases, which suggests the existence of an "optimum" optical geometry 
wherein radial gradients are minimized and total-power deposition is maximized. Such 
an optimum geometry might be achieved by increasing dye concentration and focusing 
solar radiation more toward the center to compensate for attenuation in outer regions 
of the tube. This calls for a particular concentrator design for each particular dye 
concentration and geometry, and the resultant distribution of power deposition would 
be valid only for the unique directional, spectral, and emissive properties of the 
Sun. 


Calculations show quantitatively that improvements on the postulated model are 
necessary to reach lasing threshold. Several possible improvements are identified 
and their effects are estimated from data here and in referenced literature. The 
effects are estimated conservatively and indicate that the active length of the model 
laser system must be 4.72 m (at a solar concentration of 20 000) to achieve lasing 
threshold. Since it is not possible to achieve the required solar concentration over 
the required length in practice, it is not possible to build a solar-pumped laser 
with the geometry of the model. Another geometry that provides the required length 
and concentration (with essentially unchanged laser-cavity characteristics) might be 
possible if lasing threshold actually varies with length as theory predicts. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 11, 1984 
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APPENDIX 


LABORATORY MEASUREMENTS OF DYE-LASER MODEL 

The laboratory laser model consisted of a xenon flash lamp and dye-laser tube 
housed inside an almost cylindrical reflective cavity. (Actually, the cavity cross 
section was an almost circular ellipse with the laser tube at one focus and the flash 
lamp at the other focus.) The flash lamp was driven by a capacitive network which 
simulated a transmission line and produced an approximate square wave current pulse 
of 1 -msec duration. The network was charged to voltages as large as 5 kV and was 
triggered by a pulse transformer connected in series with the lamp. 

Rhodamine-dye solution was circulated in series through a pump, connection 
tubes, laser tube, reservoir, and filter. All materials contacting the dye were 
Teflon, 1 stainless steel, and glass, except for the commercially available filter 
material that was unknown. 

The laser tube was a cylindrical Pyrex glass tube 304.8 mm long with an inside 
diameter of 4.0 mm and an outside diameter of 6.35 mm filled with rhodamine 590, 
methanol, water, and eye lo-octatetraene, an effective triplet-state quencher. The 
output mirror was 1 percent transmissive, the rear mirror was almost totally reflec- 
tive, and both had a 1 0-m radius of curvature. Dye concentration was 

1 o 

1.51 x 10 molecules per cubic millimeter. 

Energy absorbed by the dye solution was determined from volumetric expansion of 
the solution after a flash-lamp pulse and was found to be approximately 26.7 J. 

Thus, absorbed power during the 1-msec pulse was approximately 26 700 W. (Subsequent 
measurements with an improved system indicate that the absorbed power can be reduced 
by almost a factor of 2. ) 


teflon: Trademark of E. I. du Pont de Nemours & Co., Inc. 

2 Pyrex: Trademark of Corning Glass Works. 
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TABLE I 


INTENSITY VALUES 


a -1 

mm 

Intensity values F(a,R) at increments of R/R q 

0 

0.9600 

0.9600 

0.9600 

0.9600 

0.9600 

0.9599 

0.9598 


.9596 

.9593 

.9587 

.9577 

.9561 

.9530 

.9470 


.9334 

.8761 

.7133 

.6375 

.5832 

.5380 

.5019 

0*2 

0.6435 

0.6437 

0.6444 

0.6455 

0.6472 

0.6492 

0.6516 


• 6546 

.6579 

.6617 

.6657 

.6699 

.6737 

.6761 


.6735 

.6381 

.5196 

.4666 

.4296 

.3984 

.3730 

0.4 

0.4314 

0.4317 

0.4329 

0.4 349 

0.4376 

0.4412 

0.4456 


.4508 

.4569 

.4639 

.471 8 

.4804 

.4895 

.4984 


.5044 

.4861 

.4023 

.3676 

.3447 

.3256 

.3103 

0.6 

0.2891 

0.2896 

0.291 1 

0.2935 

0.2969 

0.301 3 

0.3068 


.3134 

.321 1 

.3301 

.3403 

.3518 

.3644 

.3778 


.3899 

.3842 

.3271 

.3068 

.2952 

.2860 

.2792 

o 

. 

CD 

0.1 938 

0.1 943 

0. 1 958 

0.1 984 

0. 2020 

0.2067 

0.2126 


.2197 

.2281 

.2380 

.2494 

.2624 

.2770 

.2930 


.3092 

.3127 

. 2759 

.2666 

.2642 

.2631 

.2633 

1.0 

0.1 299 

0.1 304 

0.1319 

0.1344 

0.1 379 

0.1425 

0.1483 


.1 553 

.1 638 

.1 737 

.1 854 

.1988 

.2143 

.231 7 


.2502 

.2601 

.2386 

.2380 

.2430 

.2488 

.2554 

1.2 

0.0871 

0.0875 

0.0889 

0.091 2 

0.0944 

0.0987 

0.1 040 


.1107 

.1186 

.1282 

.1394 

.1527 

.1681 

.1859 


.2056 

.2199 

.2099 

.2161 

.2269 

.2388 

.2518 

1.4 

0.0584 

0.0588 

0.0599 

0.0619 

0.0648 

0.0686 

0.0734 


.0794 

.0867 

.0954 

.1 060 

.1 185 

.1 334 

.1 510 


.1710 

.1881 

.1868 

.1982 

.2140 

.2312 

.2498 

1.6 

0.0391 

0.0395 

0.0405 

0.0422 

0.0446 

0.0479 

0.0521 


.0573 

.0638 

.071 7 

.081 3 

.0929 

.1069 

.1237 


.1 435 

.1 623 

.1 674 

.1 829 

.2028 

.2247 

.2485 

1.8 

0.0262 

0.0265 

0.0273 

0.0287 

0.0308 

0.0335 

0.0371 


.0416 

.0472 

.0542 

.0628 

.0733 

.0863 

.1 022 


.121 3 

.141 0 

.1508 

.1694 

.1928 

.2190 

.2479 

2.0 

0.0176 

0.0178 

0.0185 

0.0196 

0.0213 

0.0236 

0.0266 


.0304 

.0352 

.041 2 

.0488 

.0582 

.0701 

.0848 


.1 031 

.1 230 

.1 362 

.1 573 

.1836 

.2137 

.2476 




Figure 1 . - Geometry and variables 




Figure 2.- Intensity distribution in dye 
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Figure 3.- Absorbed-power density in dye. 


Relative absorbed-power density 


1.0 







1. Report No, 2. Government Accession No. 

NASA TP-2321 ‘ 

4. Title and Subtitle 

INTENSITY AND ABSORBED-POWER DISTRIBUTION 
IN A CYLINDRICAL SOLAR-PUMPED DYE LASER 


7. Author(s) 


3. Recipient's Catalog No. 

5. Report Date 

May 1 984 

6. Performing Organization Code 

506-55-73-01^ 

8. Performing Organization Report No. 


M. D. Williams 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

15. Supplementary Notes 


L-15694 
10. Work Unit No. 


11. Contract or Grant No. 


13. Type of Report and Period Covered 

Technical Paper 

14. Sponsoring Agency Code 


16. Abstract 

The internal intensity and absorbed-power distribution of a simplified hypothetical 
dye laser of cylindrical geometry is calculated. Total absorbed power is also 
calculated and compared with laboratory measurements of lasing-threshold energy 
deposition in a dye cell to determine the suitability of solar radiation as a pump 
source or, alternatively, what modifications, if any, are necessary to the hypo- 
thetical system for solar pumping. 


17. 


Key Words (Suggested by Author(s) ) 


18 . 


Distribution Statement 


Dye lasers 
Solar pumping 
Energy conversion 


Unclassified 


Unlimi ted 


Subject Category 36 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price 

Unclassified 

Unclassified 

16 

A0 2 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1984 



THIRD-CLASS BULK RATE 


National Aeronautics and 
Space Administration 

Washington, D.C. 

20546 

Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 



2 1 1 U ,D , 840511 S00PI6DS 

DEPT OF THE AIR FORCE 
AF SYSTEMS COMMAND 

ATTN* AFSC LIAISON OFFICE JV) S 2.2 1 
NASA LANGLEY RESEARCH CENTER 
LANGLEY AFB VA 23665 




POSTMASTER: 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 




